Specific Heat of Mg n B 2 
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Measurements of the specific heat of Mg 11 B2, from 1 to 50 K, in magnetic fields to 9 T, give the 
Debye temperature, = 1050 K, the coefficient of the normal-state electron contribution, y n = 2.6 
mj mol -1 K -2 , and a discontinuity in the zero-field specific heat of 133 mj mol -1 K _1 at T c = 
38.7 K. The estimated value of the electron-phonon coupling parameter, A = 0.62, could account for 
the observed T c only if the important phonon frequencies are unusually high relative to 0. At low 
T, there is a strongly field-dependent feature that suggests the existence of a second energy gap, 
about four times smaller than the major gap. 

PACS numbers: 74.25.Bt, 74.60.Ec, 75.40.Cx 



The superconductivity of MgB 2 , with T c = 39 K @] 
and an isotope effect ^, [5| Q consistent with the phonon- 
mediated electron pairing of the BCS theory |5), has re- 
opened the question of the maximum T c that can be pro- 
duced by that mechanism || 0]. It also raises the com- 
plementary question: What is the mechanism of the su- 
perconductivity of MgB2 itself? In this Letter, we report 
measurements of the heat capacity (C) of Mg n B2 that 
give information relevant to the latter question: Com- 
parison of the coefficient (7„) of the normal-state elec- 
tron contribution to C with band-structure calculations 
H ||, |l(J gives an estimate of the electron-phonon cou- 
pling parameter (A). The value of A suggests moderately 
strong coupling, but whether it can account for the high 
value of T c on the basis of phonon coupling alone depends 
on the values of other parameters that are not yet deter- 
mined. In the superconducting state C deviates from 
the BCS expression in a way that has no parallel among 
known superconductors, and which suggests the presence 
of a second, smaller gap in the electron density of states 
(EDOS), which has a non-BCS dependence on T. 

The Mg n B 2 sample was a powder, prepared by react- 
ing 11 B powder and Mg metal in a capped BN crucible 
at 850°C under a 50-bar argon atmosphere for 1.5 hours. 
Thermal contact to the powder was achieved by mixing it 
with a small amount of GE7031 varnish (a common low- 
T thermal contact agent with a known heat capacity) in 
a thin-walled copper cup. These extra contributions to 
the addenda limited the precision of the data. However, 
the alternate method of providing thermal contact, sin- 
tering the powder, can have adverse effects on the sam- 
ple |Il| and may in some cases account for differences 
between the results reported here and those obtained in 
other measurements. The measurements of C were made 
by a modified heat-pulse technique, 1-32 K, and by a 
continuous-heating technique, 29-50 K. Measurements 
in magnetic field (H) were made on the field-cooled (FC) 
sample after applying the field at T > 60 K. For H = 1 
T, for which the field penetration is about 50% of that 
in the normal state, measurements made after applying 
the field at 1 K were indistinguishable from the FC re- 
sults, suggesting that equilibrium flux distributions were 



attained. 

Below 2 K, there is an if-dependent hyperfine con- 
tribution to C . There are also several //-independent 
features in C, including an "upturn" below 2 K, that 
are probably associated with small amounts of impurity 
phases (or possibly with the GE varnish, although past 
experience makes that seem unlikely). Partly for that 
reason, most of the interpretation of the results is based 
on an analysis of the differences, C(H) — C(9 T), in which 
the if-independent extraneous contributions and most of 
the contributions of the addenda, including the varnish, 
cancel. C(H) — C(9 T) was calculated after the data were 
corrected for the hyperfine contributions and a small in- 
dependent part of the heat capacity of the sample holder. 
Since C is the sum of an iJ-dependent electron contri- 
bution (C e ) and an //-independent lattice contribution 
(Ci), an analysis of C(H) — C(9 T) also has the advan- 
tage that Ci cancels, leaving C e (H), the contribution of 
greater interest. In the normal state C e (H) = 7„T, in- 
dependent of H; in the mixed state C e (H) includes a 
T-proportional term, ~f(H)T, and //-dependent terms; 
in the superconducting state, C e (0) = C es 

The upper critical field (H C 2) of MgB 2 is approximately 
linear in T with H c2 (0) ~ 16 T |2). This is reasonably 
consistent with the C measurements, for which the onset 
of the transition to the mixed state, at H — H C 2(T), is 
marked by the deviations of C(H) — C(9 T) from zero 
(see Fig 1(a)). It leads to the expectation that for 9 T 
the sample would be in the normal state for T > 20 K. 
The data actually plotted in Fig. 1(a) are the differ- 
ences, [C{H) - (7(9 T)]/T, but with the scale shifted by 
7(9 T). (They are essentially point-to-point differences, 
with no smoothing of the 9-T data, which increases the 
scatter.) The values of y(H) have been determined by fit- 
ting the low-T, mixed- and superconducting-state data 
with [C(H)-C(9 T)] = [7(^-7(9 T)]T+aexp(-6/T), 
where a and b are H dependent. Consistent with the 
H = data, 7(0) was taken as zero, fixing the values 
of j(H) for all H. Figure 2 shows j(H) vs H, and the 
extrapolation to -ff C 2(0) = 16 T to obtain j n = 2.6 mj 
mol -1 K -2 . Although the extrapolation is somewhat ar- 
bitrary, the very small differences between 7(5 T), 7(7 
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T), and 7(9 T) suggest that it gives j n to within ~ 0.1 
mJ mol -1 K~ 2 . If the sample were normal at all temper- 
atures in 9 T the quantity shown in Fig. 1(a) would be 
exactly C e (H)/T. That is not the case, but the differ- 
ences are small, as shown by the small differences between 
C*(5 T), C(7 T), and C*(9 T). Quantitatively, C e (H)/T 
is underestimated by the amount j n — C e (9 T)/T. For 
T < 10 K, 7n - C e (9 T)/T « 7n - 7(9 T) « 0.08 mJ 
mob 1 K~ 2 ; for 10 < T < 20 K, where there must be 
a small broad anomaly in C(9 T), the underestimate is 
smaller and T-dependent. 

The transition to the superconducting state, shown in 
Fig. 1(c), with an entropy-conserving construction that 
gives T c = 38.7 K and AC(T C ) = 133 mJ mob 1 K"\ 
is relatively sharp, with a width ~ 2 K. In this temper- 
ature interval the sample is in the normal state for H 
= 9 T, and addition of 7 „ = 2.6 mJ mol" 1 K~ 2 to the 
quantity plotted would give the electron contribution to 
C through the transition. The effect of H in broadening 
the transition (to the mixed state), as expected for mea- 
surements on a powder with an anisotropic H C 2 ]9|, [l3| , 
is evident in Fig. 1(a). 

The thermodynamic consistency of the data, including 
in particular the very unusual T dependence of C es (0), 
can be tested by calculating the difference in entropy (S) 
between and T c for different fields. The result of such 
a test is shown in Fig. 3(a) where the entropies ob- 
tained by integrating the plotted points are compared 
with 7(9 T)T, which represents the 9 T data. At 40 
K, the entropies for all H are within ± 2% of the same 
value. The result of a second integration of the entropies 
to obtain free energy differences and the thermodynamic 
critical field (H c ) is shown in Fig. 3(b). 

Fitting the 9-T data for 20 < T < 50 K with (7(9 T) = 
7„T + C h where d = B 3 T 3 + B 5 T 5 , gave B 3 = 5.1 x 
10~ 3 mJ mob 1 K~ 4 and B 5 = 2.5 x 10~ 6 mJ mob 1 
K~ 6 . The Debye temperature (G), calculated following 
the usual convention of using the value of B 3 per g atom, 
is 1050 ± 50 K. 

Fig. 1(a) includes a comparison of the experimental 
C es with that for a BCS superconductor with j n = 2.6 
mJ mol" 1 K" 2 and T c = 38.7 K. For T > 27 K, C es is 
approximately parallel to the BCS curve; at lower T it 
rises above the BCS curve and then decreases to zero as 
aexp(— b/T), but with values of a and b very different 
from those of the BCS superconductor. It seems unlikely 
that this behavior could be understood on the basis of 
gap anisotropy similar to that known in other supercon- 
ductors. Qualitatively, it gives the appearance of a tran- 
sition to the superconducting state in two stages: the 
first, a partial transition at T c that leaves a "residual" 
7 (the extrapolation of C es to T = from above 12 K 
gives ~ 1 mJ mol -1 K~ 2 , but it should be regarded as 
an overestimate because it does not allow for the full en- 
tropy associated with the second stage); a second stage 
that is associated with a second, smaller, energy gap, 
which decreases in amplitude in the vicinity of 10 K. Tak- 
ing the exponential decrease in C es as a manifestation of 



a BCS-like transition, and comparing the parameters a 
and b with BCS theory, gives T c — 11 K and 7 = 0.74 mJ 
mob 1 K~ 2 . This interpretation of C es is, to some degree, 
understandable on the basis of theoretical considerations: 
It has been suggested that the gaps may be different on 
two parts of the Fermi surface |L4| , and, depending on 
the strength of the electron-phonon coupling between, 
and on, the two parts, the amplitude of the small gap 
can show a relatively abrupt decrease at a temperature 
well below T c . Although both gaps open at T c , that fea- 
ture in the small gap could produce the observed feature 
in C es in the 8-12 K region. The existence of two gaps on 
the Fermi surface is also consistent with scanning tunnel- 
ing spectroscopy, which has shown both a flat-bottomed 
BCS-like gap at low T, but with a small amplitude, ~ 2 
meV, corresponding to a BCS T c of ~ 13 K fL5|| , and 
a V-shaped gap with an amplitude of ~ 5.2 meV, cor- 
responding to a BCS T c of 35 K |l6|], more compatible 
with the observed T c . Despite the difficulty of explain- 
ing why and how different groups measure different gaps 
(see also |l7j and Jl8|]), it is striking how well these two 
gaps would account for C es : At low T, where the thermal 
excitations are too weak to overcome the larger gap, the 
exponential behavior of C es is consistent with the smaller 
gap; the opening of the large gap at T c explains the large 
AC(T C ). Moreover, between 20 K and T c , C es /T has a 
linear behavior (see Fig. 1(a)), which is consistent with 
a V-shaped gap. (Similar T 2 behavior has been seen in 
heavy-fermion superconductors |1£].) 

Anisotropy in H C 2 cannot explain the dramatic in- 
crease in 7 (if) at low H shown in Fig. 2. The dashed 
curve is a calculation using the effective-mass model with 
an anisotropy of 10, which is already greater than re- 
ported values (|, but 7 (if) cannot be fitted with 
any value of the anisotropy. The rapid increase in j(H) 
at low H reflects the difference between C e (0.5 T) and 
C e (0) below 10 K, and is related to the existence of the 
second gap. 

An average over the Fermi surface of the electron- 
phonon coupling parameter can be estimated by com- 
paring 7„ with band-structure calculations of the "bare" 
EDOS at the Fermi surface (N(0)) using the relation 
7„ = (l/3)7r 2 fc 2 3 A^(0)(l + A). N(0) has been reported 
as 0.68, 0.71, and 0.72 states eV" 1 unit cell" 1 (Refs. [§, 
§, and H respectively) giving A = 0.62, 0.56, and 0.53. 
Theoretically calculated values of A are 0.68 ||, and ~ 1 
0. In the Mc Millan relation [§, T c is related to A, 6, 
and the electron electron repulsion (fi*) by 

T c = (9/1.45) exp{-l. 04(1 + A)/[A - //(l + 0.62A)]}, 

(1) 

with ^i* frequently taken to be ~ 0.1 |J, [|. With 
A = 0.62, the highest of the values derived from 7„ and 
A(0), and 6 = 1050 K Eq. (1) gives T c = 22 K, too low 
by almost a factor two. With these values of A and 0, 
T c = 39 K would require fx* — 0.030, an unusually low 
value. However, in this expression O represents a rela- 
tively crude estimate of the phonon frequencies that are 
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important in the electron pairing. As derived from the 
coefficient of the T 3 term in C/ , 8 is really a measure of 
the frequencies of the low-frequency acoustic phonons, 
which may not be particularly relevant to the pairing of 
the electrons. In the Debye model, 8 is also the cut-off 
frequency but real phonon spectra often extend to signifi- 
cantly higher frequencies. In the Allen and Dynes version 
of the theory 0, 6/1.45 is replaced by w; og /1.20, where 
ujiog is a moment of the phonon frequencies in which they 
are weighted by the electron-phonon matrix elements. 
More detailed calculations that take into account relevant 
features of the phonon spectrum and electron-phonon 
scattering may give a value of the pre-exponential factor 
in Eq. (1) that accounts for the observed T c with a phys- 
ically plausible value of /x* , but until they are available, 
the question of whether the electron pairing in MgB2 is 
phonon mediated would seem to remain open. 

The parameters (H c (0)) 2 /j n T* and AC[T C )/ 7 „T C 
measure the strength of the electron pairing P0| . In the 
BCS, weak-coupling limit, their values are 5.95 and 1.43, 
respectively. There are a number of "strong-coupled" 
superconductors for which these parameters are greater 
than the BCS values, but relatively few for which they are 
smaller |2^]. For Mg 11 B2 they are unusually small, 5.46 
and 1.32, which, in contrast with the value of A, suggests 
extreme weak coupling. This indication of weak coupling, 
based on bulk properties, is in qualitative agreement with 
tunneling results that may have been influenced by sur- 
face effects Q. 

In general, there are more differences than similar- 
ities among the specific heat measurements on MgB2 
^ EHl |2^, |2^, However, in several important re- 

spects, including the low— T behavior of C es , the results 
of the Geneva group EJ are qualitatively similar to those 
reported here, and the differences that do occur are read- 
ily understood in terms of sample dependence. A com- 
parison with their results, which were obtained by differ- 
ent experimental techniques on sintered sample of com- 
mercial material, attests the qualitative validity of the 
major features reported both here and in Ref. The 
value of 7„, 2.7 mJ mol -1 K~ 2 , and the test of ther- 
modynamic consistency in Ref. j2^] are similar and of 
comparable accuracy to those reported here, but limited 
by paramagnetic Fe impurities in Ref. pl| and the pre- 
cision of the data here. Although the superconducting- 
state entropies are essentially identical at T c , the specific 
heat anomaly at T c in Ref. 21 1 is broader, leading to a 
substantial underestimate of AC(T C ). Most importantly, 
the T-dependence of C es in the 4-15 K range, which 
shows the presence of a second gap (or extreme an un- 
usual anisotropy), is essentially the same in both cases. 
However, for T < 4 K, C es was obscured by the contribu- 
tion of paramagnetic impurities and the limiting T — » 
dependence was taken as approximatively T 2 |2l| ] , which 
would be expected for line nodes, rather than the expo- 
nential dependence reported here. 

Several other values of AC(T C ), all for sintered samples 
and all lower than that reported here, have been given 



in other reports [g, [2^, E3L |24j . Values of j n that range 
from 1.1 to 5.5 mJ mol~MC , based on different [^2[ |2J] 
or unspecified [|| analysis of experimental data have also 
been reported, but the very similar values reported here 
and in Ref. M, are supported by the thermodynamic 
consistency of the data 
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FIG. 1: (a) [C(H) - C(9 T)]/T. In (a) and (b) the scale has 
been shifted by 7(9 T) to give an approximation to C e (H)/T 
(see text). In (a) the dashed curve is a polynomial extrapo- 
lation of the 12-20 K, H = data to T = 0; the horizontal 
line represents 7(9 T). In (b) the low-T 5- and 7-T data 
are shown on an expanded scale with solid curves represent- 
ing fits described in the text; the horizontal solid and dashed 
lines represent 7(9 T) and y„, respectively. In (c) the solid 
lines represent an entropy conserving construction. The error 
bars are ± 0.1% C to tai/T. 



FIG. 2: 7 as a function of H. The solid curve is a guide to 
the eye and an extrapolation to H C 2(0). The dashed curve is 
a "fit" with an H C 2 anisotropy of 10. 



FIG. 3: (a) Entropies as functions of T for different H, with 
H increasing from the lowest to highest curve, (b) Thermo- 
dynamic critical field, compared with a BCS curve for the 
derived values j n and T c . 



